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CN ■ In this paper, we propose an approximate projected consensus algorithm for a network 

to cooperatively compute the intersection of convex sets. Instead of assuming the exact 
convex projection proposed in the literature, we allow each node to compute an approxi- 
mate projection and communicate it to its neighbors. The communication graph is directed 
\^ , and time-varying. Nodes update their states by weighted averaging. Projection accuracy 

04 ' conditions are presented for the considered algorithm. They indicate how much projection 

accuracy is required to ensure global consensus to a point in the intersection set when the 
, communication graph is uniformly jointly strongly connected. We show that 7r/4 is a crit- 

' ical angle error of the projection approximation to ensure a bounded state. A numerical 

^ . example indicates that this approximate projected consensus algorithm may achieve better 

O I performance than the exact projected consensus algorithm in some cases. 
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^ ' 1 Introduction 
in 
O 

. In recent years, dynamics and control on large-scale networks have drawn increasing research 

attention in different areas including engineering, computer science, and social science. Co- 
operative control of a group of autonomous agents fully employs local information exchange 
^ \ and distributed protocol design to accomplish collective tasks such as consensus, formation, 

^ ; and aggregation [ZlElliaiEllIlllMllISlinillllllS]- Moreover, in parallel computation, load- 

balance problems require realtime balance of the load from different computing resources [9l |T0] . 
Additionally, a central problem of opinion dynamics in social networks is how the agreement 
is achieved via individual belief exchange processes [131 114j . A fundamental question in these 
problems is, how consensus can be guaranteed based on local information exchange, time- varying 
node interconnection and limited knowledge of the global objective. 
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Various distributed optimization problems arise for consensus with particular optimization 
purpose in practice. Minimizing a sum of convex functions, where each component is known 
only to a particular node, has attracted much attention, due to its simple formulation and 
wide appHcations [SaiMlEIllMlEZllSaiMliniEaiMlEaEaiSl]. The key idea is that 
properly designed distributed control protocols or computation algorithms can lead to a collec- 
tive optimization, based on simple exchanged information and individual optimum observation. 
Subgradient-based incremental methods were established via deterministic or randomized itera- 
tion, where each node is assumed to be able to compute a local subgradient value of its objective 
function [20l EU [26l [22l [25l [2l] . Non-subgradient-based methods also showed up in the liter- 
ature. For instance, a non-gradient-based algorithm was proposed, where each node starts at 
its own optimal solution and updates using a pairwise equalizing protocol [291 [30], and later an 
augmented Lagrangian method was introduced in [33]. 

In particular, if the optimal solution set of its own objective can be obtained for each node, 
the considered optimization problem is then converted to a set intersection computation problem 
when we additionally assume there is a nonempty intersection among all solution sets [32l [311 [27] . 
In fact, convex intersection computation problem is a classical problem in the optimization study 
[351 |36l [37] . The so-called "alternating projection algorithm" was a standard centralized solution, 
where projection is carried out alternatively onto each set [35l [361 [37] . Then the "projected con- 
sensus algorithm" was presented as a decentralized version of alternating projection algorithm, 
where each node alternatively projects onto its own set and averages with its neighbors, and 
comprehensive convergence analysis was given for this projected algorithm under time-varying 
directed interconnections in [27J. Following this work, a flip-coin algorithm was introduced when 
each node randomly chooses projection or averaging by Bernoulli processes, and almost sure con- 
vergence was shown for the system to reach an optimal consensus in [32]. A dynamical system 
solution was given in [31], where the network reaches a global optimal consensus by a simple 
continuous-time control. In all these algorithms, each node needs to know the exact convex 
projection of its current state onto its objective set [32l [3T1 [27]. 

However, in practice, the exact convex projection is usually hard to compute due to the 
common environmental noise and computation inaccuracy. In this paper, we therefore propose 
an approximate projected consensus algorithm to solve the convex intersection computation 
problem. Instead of assuming the exact convex projection, we allow each node to just compute 
an approximate projection point which locates in the intersection of the convex cone generated 
by the current state and all directions with the exact projection direction less than some angle 
and the half-space containing the current state with its boundary being a supporting hyperplane 
to its own set at its exact projection point onto its set. The communication graph is supposed 
to be directed and time-varying. With uniformly jointly strongly connected conditions, we show 
that the whole network can achieve a global consensus within the intersection of all convex sets 
when sufficient projection accuracy can be guaranteed. For a special approximate projection 
case when the nodes can get the exact direction of the projection, a necessary and sufficient 
condition is given on how much projection accuracy is critical to ensure a global intersection 
computation. A numerical example is also given, and surprisingly, the approximate projected 
consensus algorithm sometimes achieves better performance for convergence than the exact 
projected consensus algorithm. 
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The paper is organized as follows. Section 2 gives some basic concepts on graph theory and 
convex analysis. Section 3 introduces the network model and formulates the problem of inter- 
est. Section 4 presents the main results and convergence analysis for the considered algorithm. 
Section 5 focuses on the discussion on the critical angle error of the approximate projection. 
Section 6 gives a numerical example and finally, Section 7 shows some concluding remarks. 

2 Preliminaries 

In this section, we introduce preliminary knowledge on graph theory [5] and convex analysis [1]. 

A directed graph (digraph) G = (y,£,A) consists of node set V = {1,2, arc set 

£^ C V X V and an adjacency matrix A = [aij]nxn with nonnegative adjacency elements aij. 
The element aij of matrix A associated with arc is positive if and only if £ £■ Mi 
denotes the set of neighbors of node i, that is, Mi = {j S V\{i,j) S £}. In this paper, we assume 
G £ for all i. A path from i to j in digraph ^ is a sequence (io,«i), {11,12), (V-ii V) °f 
arcs with iq = i and ip = j. Q is said to be strongly connected if there exists a path from i to j 
for each pair of nodes i, j G V. 

A function /(•) : BJ^ ^ R is said to be convex if /(Ax + (1 - \)y) < \f{x) + (1 - \)f{y) 
for any x,y G BJ^ and < A < 1. A function / is said to be concave if — / is convex. A set 
K C i?™ is said to be convex if Ax + (1 — \)y G K for any x,y € K and < A < 1 and is 
said to be a convex cone if Aix + X2y G K for any x,y £ K and Ai,A2 > 0. For a closed 
convex set K in i?™, we can associate to any x G a unique element Pk{x) G K satisfying 
\x — Pk{x)\ = miy^K \x — y|, which is denoted as \x\k, where | • | denotes the Euclidean norm 
and Pk is the projection operator onto K. 

For a closed convex set \i x ^ K, then by the supporting hyperplane theorem, there 
is a supporting hyperplane to K at Pk{x). The angle between vectors a and h is denoted as 
Ang(a,6) G [0, tt] for which cosAng(a, 5) = (a, 6)/(|a| where (a, 6) denotes the Euclidean 
inner product of vectors a and h. 

We cite a lemma for the following analysis (see Example 3.16 in [3], pp. 88). 

Lemma 2.1 f{z) = \z\k is a convex function, where K is a closed convex set in R^. 

The following properties hold for the projection operator Pk- Here (i) is the standard non- 
expansiveness property for convex projection; (ii) comes from exercise 1.2 (c) in [2] (pp. 23) and 
(iii) is a special case of proposition 1.3 in \T\ (pp. 24). 

Lemma 2.2 Let K be a closed convex set in i?™. Then 

(i) iPxix) - PK{y)\ <\x-y\ for any x andy; 

(ii) \\x\k — \y\K\ ^ \x — y\ for any X and y; 

(Hi) Pk{Xx + (1 — X)Pk{x)) = Pk{x) for any x and < A < 1. 

The next lemma can be found in [32] . 
Lemma 2.3 Let K and Kq <^ K be two closed convex sets. We have 

\Pk{x)\k^ + \x\k < \x\k^^ for anyx. 
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3 Problem Formulation 



In this section, we introduce the intersection computation problem and the approximate pro- 
jected consensus algorithm. 

Consider a network consisting of n agents with node set V = {l,2,...,n}. Each node i is 
associated with a set Xi C R"^ and set Xi is known only by node i. The intersection of all these 
sets is nonempty, i.e., HILi -^i ®- ^et us denote Xq = HILi -^i- "^^^ target of the network is 
to find a point in Xq in a distributed way. For Xi, i = 1, n, we use the following assumption: 

Al (Convexity) Xi,i = 1, ...,n, are closed convex sets. 

Remark 3.1 The intersection computation problem can be equivalently converted into the fol- 
lowing distributed optimization problem: the objective of this group of n agents is not only to 
achieve a consensus, but also to cooperatively solve 

n 

1=1 

where fi : i?™ R is the convex cost function of agent i and can be known only by agent i. 
Here Xi = {y\fi{y) = miiix^R^ fi{x)},l < i < n are closed convex sets, which are assumed to 
be nonempty and have a nonempty intersection. 

3.1 Communication Graphs 

The communication over the network is modeled as a sequence of directed graphs, Qk = 
{V , £ (k) , A{k)) , k > 0. We say node j is a neighbor of node i at time k if there is an arc 
{i,j) G £{k), where aij{k) represents its weight. Let Mi{k) denote the set of neighbors of agent 
i at time k. We introduce an assumption on the weights [261 132] . 

A2 (Weights Rule) (i) Ylij(^Ni(k) '^iji^) — 1 ^ ^'^'^ ^■ 

(ii) There exists a constant < r/ < 1 such that aij{k) > r] for all i, k and j G Mi{k). 

For the connectivity of the communication graphs, we introduce the following assumption 
[311 [27]- 

A3 (Connectivity) The communication graph is uniformly jointly strongly connected (UJSC), 
i.e., there exists a positive integer T such that Q{[k,k + T)) is strongly connected for A; > 0, 
where G{[k, k + T)) denotes the union graph with node set V and arc set Ufc<s<fc+T'^(^)- 

3.2 Approximate Projection 

Projection methods have been widely used to solve various problems, including projected con- 
sensus [27], the convex intersection computation |36l [37] and distributed computation [3]. In 
the most literature, the projection point Pk{z) of z onto closed convex set K is required to 
achieve desired convergence, but in practice it is hard to be obtained and often is computed 
approximately. Here is the definition of approximate projection. 
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Definition 3.1 Suppose K C i?™- is a closed convex set and < 9 < Tr/2. Define 

Ck{v, 9) = V + {z\ {z, Pk{v) — v) > \z\\v\k cos 0}; 
U+{v) = {z\ {v-Pk{v),z) > {v-Pk{v),Pk{v))}. 

as a convex cone and a half-space, respectively. The approximate projection of point v 

onto K with approximate angle 9 is defined as the following set: 



Ck{v,9)^U+{v), ifv^K 
{v}, if V €z K 



(1) 



In fact, Ck{v, 9) — v is a. convex cone generated by all vectors having angle with Pk{v) — v 
less than 9 and H^(f) is the half-space containing point v with 

YLk{v) := {z\ {v-Pk{v),z) = {v-Pk{v),Pk{v))] 
being a supporting hyperplane to K at Pk{v). 




Figure 1: The set marked by the shaded area is the approximate projection of v onto K. 



Definition 3.2 The supporting approximate projection ^f^(v,9) of point v onto K with ap- 
proximate angle 9 is defined as 

^ ' \{v}, if vGK 

Remark 3.1 Exact projection may be difficult to obtain in practice, due to possible data quan- 
tization and limitation of the estimate. This is why we introduce approximate projection, as a 
relaxed estimate to the real projection value. 

According to Definition 13.21 for any y £ we can associate y with y G ^f^{v,9) such 

that 

y = (1 - /3)t; + Py for some < /3 < 1. (2) 
Moreover, it is easy to see that \i y ^ v, y satisfying ([2]) is unique. 
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3.3 Distributed Iterative Algorithm 

To solve the intersection computation problem, we propose the following approximate projected 
consensus algorithm: 

Xi{k + l)= a,j{k)P^^{k), i = l,...,n, (3) 

where P^{k) G (xj(/c), 6*^) for all i and k. 

According to the definition of supporting approximate projection, there exist < Oi^k ^ 1 
and Pf^'ik) e ^'£^{xi{k),ek) such that 

P^ik) = (1 - ai,k)xi{k) + ai^kPrik). (4) 

Combining with 1^ and (HD, we have 

Xiik + 1)= a,,{k)(^il-aj,k)xjik) + aj,kPrik)), (5) 

where if x^{k) Xi, Pl^{k) G UxM{k)) and Ang{Pf^{k) - x,{k), PxMik)) - Xi{k)) < Ok- 




Figure 2: The approximate projected consensus algorithm. 

Remark 3.2 Note that the "projected consensus algorithm" presented in \27^ is a special case 
with Oi^k = 1 o,nd 9^ = of the approximate projected consensus algorithm discussed in this 
paper. 

We illustrate the iteration process of the algorithm in Figure 2. For the approximate angle 
9k, we use the following assumption. 
A4 < 6ife < 6** < 7r/2 for all k. 

In this paper, we are interested in whether an optimal consensus can be achieved or not, as 
defined in the following definition. 
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Definition 3.3 A global optimal consensus is achieved for the approximate projected consensus 
algorithm if, for any initial condition x{0) £ i?"™", there exists x* £ Xq such that 

lim Xi{k) = X* , i = I, n. 

Note that a global optimal consensus x* necessarily belongs to the intersection set Xq. Hence, 
a possible algorithm to find a point in Xq is to employ the approximate projected consensus 
algorithm ^ . We next discuss the convergence of this algorithm. 

4 Main Results and Convergence Analysis 

In this section, we present the main results and convergence analysis. 
Denote = mini<i<„ Oj^fc and = maxi<j<„ Oi^k, k>0. 

Theorem 4.1 Suppose A1-A4 hold. Global optimal consensus is achieved for the approximate 
projected consensus algorithm ^/ X^fcLo '^fc ~ °° ^^'^ Yl'k^o'^k^k < oo. 

To investigate the necessity of the divergence condition in Theorem 14. we impose another 
assumption on the boundedness of the n sets Xi,i = 1, n. 
A5 (Bounded Sets) Xi,i = l,...,n, are bounded sets. 

Then the following conclusion holds showing a necessary projection accuracy condition. 

Theorem 4.2 Suppose A1-A5 hold, 6^ = and < 1 for all k. Global optimal con- 
sensus is achieved for the approximate projected consensus algorithm only if X^^Lo '^fc ~ 
In fact, if '^'kLoOik ^ ^'^^ then for initial condition Xj(0) = z*,i = l,...,n with \z*\xo > 
T,kLo(^kd*/U.'k=oi'^ - "fe)' ^^ere is y* = y*{z*) Xq such that \\m.k^^Xi{k) = y* for all i, 
where d* = sup^^^^^eUILi ^« '^i ~ '^2|- 

Remark 4.1 Suppose A1-A5 hold, 6k = and there exists a sequence {ofel^o '^^^^ c^fc < 1 
for all k such that ai^k = ctj,k = c^fc for all i,j and k. Then from Theorems 1 and 2, we have 
that global optimal consensus is achieved for the approximate projected consensus algorithm if 
and only if Ofc = oo. 

Remark 4.2 Gompared to the convergence results given in \27^ , Theorems and \4-S\ do not 
require the doubly stochastic assumption on the weights aij{k) (Ylj=i'^iji^) — Ylj=i^jii^) ~ ^ 
for all i,k). This is important because double stochasticity is hard to guarantee for the arc 
weights in a distributed way, especially when the communication is directed. 

Moreover, the connectivity assumption in \21^ requires that Q{[k,k + T)) is a fixed graph 
for sufficiently large k, which is more restrictive than our UJSC assumption. However, the 
assumption in ^27^ can be relaxed to UJSG graphs, as indicated in a comment by the authors. 

4.1 Lemmas 

We establish several useful lemmas in this subsection. 
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Let {xi{k)}'^^Q be the states of node i generated by ([3]), i = l,...,n. Denote \x{k)\xo 



i\xiik)\xo---\xnik)\xoV , yik) = iyi{k)...yn{k))^ with yi{k) = \xi{k)\xo- ^J\xi{k)W^ - \xi{k)W,. 
Denoting Dj. = diagjai fc, a2,fc) •••) On.fc}) the following lemma holds. 

Lemma 4.1 Suppose Al holds. For all k > s, we have 

\x{k + l)\xo < A{k)\x{k)\xo - A{k)Dky{k) +tanekA{k)Dk\x{k)\xo. (6) 
Proof. According to Lemma |2. 11 ([5]) implies 

\xi{k + l)\xo< Yl '^^jik){{^-a^,k)\xj{k)\xo+a^,k\P^ik)\xo)■ (7) 

By Lemma 12.21 (ii) , we have 

< \pr{k) - Px,{x,{k))\ + \Px,{xj{k))\x,. (s) 

The definition of Pj"'{k) ensures that 

|Pf(fc) - Px^{Xj{k))\ < iBJiek\xj{k)\xy (9) 

Moreover, it follows from Lemma 12.31 that for any j G V, 

\PxM(k))\x, < ^J\xJm%,-\^,m%■ (10) 

It follows from dZ]), ([HI), ([9]), and the relation \xj{k)\xj < \xj{k)\xa that 
\xi{k + l)\xo< Yl aij(A;)((l -aj-fc)kj(^)Uo + aj,k^J\xjik)\j^^ - \xj{k)\j^^^ 

+ tan6'fc ^ aij{k)aj^k\xj{k.)\xo (11) 
= Yl «»i(^)(l^i(^)Uo - OLj,k{\xj{k)\xo - yj\xj{k)\\^ - \xj{k)\\^)^ 

+ tan6'fc ^ ajj(A;)aj,fc|xj(A;)|xo- (12) 

Then the conclusion follows. □ 

It is easy to find that tan6l < {iaiie*/e*)e for < 6^ < 9*. Thus, if Y.kLo^k^k < oo, then 
Efclo"fe tan^fc < oo. 

Lemma 4.2 Suppose Al and A4 /loW. If '}^'k=o^t^k < c«, {a^i(^)}fc^o ^'^ bounded for all i. 
Proof. Lemma 1 (b) in [27] states that for closed convex K and any x G K^, 

\Pk{x) - yf <\x- yf - |x||- for all y e K. (13) 
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By taking K = Xj and z G Xq C Xj , ([T3|) leads to 

\Px,{x,{k)) -z\< ^\x,{k)-z\^-\x,{krx^. (14) 

By considering \xi{k + 1) — z\ instead of \xi{k + l)|xo) following similar procedures with ([7]), 
dH]), ([9]), and substituting (fTO]l with (fn|) . we can show that for any i £V, 



\xi{k + 1) — z\ < aij{k)\xj{k) — z\ + tan 9k aij{k)aj^k\xj{k) — z\ 

- ^ aij{k)aj^k(\xjik) - z\ - ^\xj{k) - z\^ - (15) 

By dropping the non positive term on the righthand in (jlSp . we have 

max \xi{k + 1) — 2;| < (1 + at tan^^) max \xi{k) — z\. (16) 

l<i<n l<«<n 

Therefore, for all k, we have 

k 

max + 1) — z| < IT fl + tan 0/) max |xj(0) — zl 

l<j<n ' ±1 ^ ' / l<j<n 

— /=o — 

< gEto^+tan^'i max |x,(0) - z| 

l<i<n 

< eE,So"^tan^i max |x,(0) - z\, (17) 

l<i<n 

where the second inequality follows from 1 -\- z < for z > 0. Then the conclusion follows. □ 
The next lemma is a special case of various random versions, for example, see Lemma 11 in 
[6] (pp. 50). Here we give a simple proof. 

Lemma 4.3 Let {afclfc^g '^'^^ {^fclfc^o non-negative sequences with Yl^=o^k < oo. Suppose 

flfc+i < fflfc + &A: for all k 

Then lim^._^oo Ofc is a finite number. 

Proof. Define a new sequence {ck}'^i with = — Yl^=o 'which is bounded since {ofcjfc^i 
and {Y^^Zq bi}kLi bounded. Moreover, c^+i < Cfc for k > 1 implies that lim^.^oo Cfc is a finite 
number. The conclusion follows from limj^^oo CLk = li™fc->oo Cfc + X]£o '-' 

Lemma 4.4 Suppose Al and A4 /loW. If Y^^=o^k^k < oo, i/ie following limit exists 

t} := lim max \xi{k)\xo- 

fc— )-oo l<i<ra 

Proof. Take z G Xq. Based on (fT^ and PT|) . we have 

max |xj(/c + l)|xo < max |xj(A;)|xo + tan6'fce^'=o"i^ max \xi{0) - z\. 

l<i<n l<i<n l<i<n 

The conclusion follows from the last inequality and Lemma 14.31 □ 
Denote 

??+ = limsup |xi(/c)|xo, 'Hi' = liminf |xi(A;)|xo, i G V. 

k~>oo k^oo 

Obviously, < rj^^ < rj^ < for all i. 
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Lemma 4.5 Suppose A1-A4 hold. X^^Lo "^fc ^fc ^ '^'^'^ t/iere exists some agent io ^ V 
such that T]^^ < I?, then ■!? = 0. 



Proof: Motivated by the idea of Lemma 4.3 in [3T], we prove this lemma by contradiction. 
Denote 

Qi = {ni + ^,^)/2, i G V. 

Since t?^^ < t?, there exists an increasing sequence {ki}'^^ such that (/;:/) |xo ^ f?io < ^ foi" 
/ > 0. Moreover, for any e > 0, there exist Kq = Kq{£) such that \xi{k)\xQ < + e and 
'^0 Yl'k=Ko '^t^k < e for /c > Ko and ah z, where 

do = (tanr/r) sup \xi{k)\xo, (18) 

l<i<n,fe>0 

which is finite by Lemma 14.21 Without loss of generality, we assume ko > Kq. 
Based on inequality (fT2]) . we have 

kio(^o + l)Uo < X] aiojiko)\xj{ko)\xo + aioio{ko)\xio{ko)\xo + doa'^^Oko- 

jGMo(*:o)Vo 

Therefore, \xig{ko + l)|xo < (1 - v)i'^ + e) + ??e'io + doa^^Ok^^ and then 

|xio(A:o + 2)|x„ < (1 - + +r/[(l - ?7)(i? + e) + ^0^^%,] + rjgig + doa^Q+i^'feo+i 
< (1 - r]^)ii} + £)+ r]^g^, + do ^^t^k- 

k=ko 

Similarly, we can show by induction that for r > 1, 

fco+'f— 1 

\xio{ko + r)\xo < i'i^ - + +v''Qio+do ^ atOk- 

k=ko 

Since the communication graph is UJSC, there exist agent ii 7^ io and time k^ G [kQ, k^ + T) 
such that {ii,io) G ^ik}^). As the above estimate for |xio(/co + t)\xo with Xi„(A;o)) by considering 
|xj^(/cQ + r)|xo with (/cq)|xo) we can show similarly that for r > 1, 

kl+T~l 

\x^M + r)\xn<{l-rf){^ + e)+v"\xi,{kl)\xn+do Yl '^t^k 



kl+T—l 



k=k(} 

Repeating the previous procedure on intervals [ko +pT, ko + {p+ l)T), 1 < p < n — 2, we can 
obtain nodes {12,13, ■■■,in-i} such that {ij,0 < j < n — 1} = V and 



max \xi{ko + f)\xo < (1 - V^){^ + + Qio + ^^o V a+9k 

[<i<n — ' 

< (l-r/^)(i9 + e)+r/^^io+e. 
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where T = (n — 1)T. Moreover, we can make similar analysis for /ci, /c2, ... and obtain that for 
/ > 0, 

max \xi{ki + f )|xo < (1 - V^){'& + + Qio + e> 

l<j<n 

which yields a contradiction since (1 — rj^){'d + e) + rj^ Qi^ + e < •§ provided that e is sufficiently 
small. □ 
We introduce transition matrices 

$(A;, s) = A{k) ■ ■ ■ A{s + l)A[s) for all k and s with k> s. 

Recall that r/ and T are defined in A2 and A3, respectively and T = (n — 1)T. The next lemma 
generalizes Lemma 2 in [26] on the lower bound of the entries of the transition matrices. 

Lemma 4.6 Suppose A2 and A3 hold. Then ^{k, s)ij > "rf- for all s and k > s + T — I. 

Proof: By Lemma 2 in [26], $(s + T — 1, s)ij > r/'^ for all i,j and s > 0. Moreover, according to 
A2 (i), ^i^i A{k)ii = 1 and then X^ILi ^i^i ^ + = 1 foi" ^) ^ Thus, for all i,j and 

/c > s + T- 1, 



$(fe,s)ij = (<^{k,s + f)^{s + f -l,s 



> V^>(A;,s + r)., min $(s + r-l,s) 

1=1 
f 

= V ■ 

The conclusion follows. □ 
Lemma 4.7 



n — ' V " V V n 

1=1 



2 



where vq > Vi > for all i. 



Proof: The conclusion follows from that f{z) = \/ (? — z'^ with domain [— c, c] is a concave 
function, where c > 0. □ 
Consider the following consensus model with noise Wj, 

Zi{k^\)= ^ hij{k)zj{k)^Wi{k), i = \,...,n, (19) 

where the weights bij{k),i,j £ V,k > satisfy A2. The consensus is said to be achieved for 
system ()19p if for any initial conditions, limfc_;,oo \ zi{k) — Zj{k)\ = for all 1 < i, j < n. The next 
lemma can be obtained from Theorem 1 in 



Lemma 4.8 If the communication graph of system U9\) is UJSC with \mii,^cc,Wi{k) = for all 
i, then the consensus is achieved for system i flPj) . 
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4.2 Proofs 

In this subsection, we present the proofs of Theorems 14.11 and 14.21 



4.2.1 Proof of Theorem 

Rewrite ([5]) as 

Xi{k + 1)= ^ aij{k)xj{k) + ^ aij{k)aj^k{Px,ixj{k)) - Xj{k)) 

+ aij{k)a,,k{Pr{k)-Px,{x,{k))). (20) 



Based on ([9]), the second and the third term in (|20p are not greater than 

max aik\xi{k)\xi +attan6k max \xi{k)\xi- (21) 

l<i<n l<i<n 

Note that = leads to hmfc_j.oo niaxi<j<„ \xi{k)\xi < hmfc^.oo iiiaxi<j<„ \xi{k)\xo = and 
then the term in ()2ip tends to zero as /c — t- cxd. Therefore, by applying Lemma 14.81 for ()20p . we 
have that if ■!? = 0, then the consensus is achieved. 

Moreover, we claim that if = and the consensus is achieved, then all agents will converge 
to a point in Xq. Since {xi{k)}'^^Q,i = 1, ...,n are bounded by Lemma Iir2] and the consensus is 
achieved, there is x* £ Xq and a subsequence {ki}f^^ such that limi^^ Xi{ki) = x* for all i. By 
similar analysis with (|17p . we have 

max \xi{k) - x*\ < 6^"=*; max |xi(/ci) - x*| < e^r=o °? t'^'^ max |xi(A;i) - x*| 

l<i<n l<i<n l<i<n 

for k > ki, which implies lim^,_).oo 3;j(A;) = x* for all i. 

If there exists some agent io such that r/^^ < 'd, then by Lemma |4.5| = 0. Therefore, we 
only need to prove 

= when rj^ = 'rj~ = "d for all i, 

which shall be proven by contradiction. If i? > 0, then for any e > 0, there exist Ki = Ki{e) 
such that \xi{k)\xo < '& + £ and doa^Ok < e for k > Ki and all i. We complete the proof by the 
following two steps. 

(i) Suppose -qf = t]^ = "d for all i. The consensus is achieved: lim^^oo \xi{k) — Xj{k)\ = for 
all i, j. 

Denote 

q = limsupaj,A:ki(A;)|x,, i G V. 

We prove q = for all i by contradiction. If there exists some agent io such that 
> 0, then there is an increasing time subsequence {ki}"^^ with ki > Ki such that 



12 



ih)\x, > cQq for all / and some < c < 1. Therefore, by ([11]) 

\xio{ki + l)\xo < «ioio(^)((l - a^^M)\xioik)\xo + OiioM\/\^ioik)\xo " kio(^i)lx,„ 

+ aiojiki)\xjiki)\xo+doateki 



j&Afig{ki)\io 



< r/((l - a,„,fcj(^? + e) + ^ a^.^i^ + e)^ - c\l) + (l - r,){^ + e) + e 



where do is the one in (|18|) . which yields a contradiction since the right hand side of 
is less than ■!? for sufficiently small e and sufficiently large I. 

Thus, limfc^oo 0(i^k\xi{k)\x, = for ah i. Moreover, since X^fclo ^k^k < oo, limfc^oo ct^ tan Ok 
< (tan 0*/^*) limfc_^c>o ^fc = 0. The two preceding conclusions and the boundedness of 
{xi(A;)}^Q imply that the term in ()2ip tends to zero and then the consensus is achieved 
by applying Lemma 14.81 for ()20p again. 

(ii) Suppose rjf = rj^ = ^ for all i. All agents converge to the nonempty intersection set: 
limfc^oo \xi{k)\xo = for all i. 
Denote 

n 

5 = limiuiS^ \xi{k)\xi- 

k^oo — ' 

4 = 1 

We prove that (5 = by contradiction. Otherwise, suppose 5 > 0. 
By ([S]), we obtain for /c > s, 

fc fc 



l=s l=s 
k-f+1 

= ^k,s)\x{s)\xo- Y Hk,l)Diy{l) 

l=s 

k k 

- Yl Hk,l)Diy{l) + doY,o^lGi^ (23) 



l=k-T+2 



l=s 



where T = [n — 1)T. By dropping the third term (nonpositive) on the right-hand side in 
([23]) . we obtain 



k-T+l 



\x{k + l)\x,<^{k,s)\x{s)\x,- Y Hk,l)Diy{l) + doY»t^i- (24) 

l=s l=s 

For £ = 6"^/ {An'^'d + 25), there exists sufficiently large K2 such that J2^=i \xi{k)\xi > S — e 



13 



and 'd — e < \xi{k)\xQ + £ for k > K2. For k > K2, we have 



1=1 i=l 



< n 



^ 2 
In ' 



i=l 

where the second inequahty follows from Lemma 14.71 and then 



1=1 



:=C>0. 



Namely, X^iLi^ilO — C for ^ > -^^2- Combining the preceding inequality with Lemma [ 
we have that every component of ^{k,l)Diy{l) is not less than rf^Qa^ for all K2 < / < 
k — T + 1 and all k > K2 + T — 1. Then by ()24p with taking s = K2, we obtain 



k-T+l k 

\xik + l)\xo<Hk,K2)\x{K2)\xo-r]^C E cc^^ + do api, (25) 

l=K2 l=K2 

where 1 is the vector of all ones. Observing that Y1'^K2 ~ Y^'iLk2 '^t^i ^ 
noticing lim,t^oo |3^(^)|xo = '^'^^ contradiction arises by taking the limit as A; — )• 00 in 

(EH). 



Therefore, 5 = 0, that is, there is a subsequence {ki}'^^^ such that limi_j.oo SILi \xi{ki)\xi = 
0. Since the consensus is achieved by what we have proven in the first step (i), we have 



lim \xi{ki)\x, = for ah j £ V, 

which implies i? = lim^^oo maxi<j<„ \xi(ki)\xQ = 0. 

This completes the proof. □ 

Remark 4.3 // the adjacency matrices are double- stochastic, Y1T=q '^k~'^ '^'"'^ XlfcLo '^t^k < 
00, then the convergence analysis for the optimal consensus can also he performed by the following 
procedures. 

By summing the two sides in over i = 1, ...,n, we obtain 



Y\xi{k + l)\xo < {1 + tanek)Y\xi{k)\xo -Yai^k{\xi{k)\xo - \]\xi{k)W^ - \xi{k)W^^ 

1=1 i=l j=l 

(26) 
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Summing the two sides in i26\) over k >0 and rearranging the terms, we have 

oo n 
fc=0 i=l 

oo n 

<^^ai^k{\xi{k)\xo - \]\xi{k)W^ - \xi{k)\ 



2 

X, 

k=0 i=l 

n oo n 

<^\xM\xo+^cilisii9k^\xi{k)\x^ < oo. (27) 

i=l k=0 i=l 



The assumption YlT=o '^fc ~ °° "^'^ (E2^ imply that 



i=l 

and then 



limmfY \xi{k)\x, = 0. (28) 

From ^27\) we also have 



fc— >oo 

i=l 



^lim ^a,,fc(|x,(/c)|xo - ^J\x^{k)\j,^ - |xi(/c)|2,J = (29) 

i=l 

and then for all i £V, 

limsupai^k\xiik)\xi = (30) 

k—^oo 

since if there is io £ V and {ki}f2^^ such that ajo,fcj^«o(^0lxio ^ ^ /^'^ '^^^ ^ ^^'^ some e > 0, t/ien 
for all I, 



£2 >0, 



which contradicts with i flP|) . i(;/iai /laue 6een proved in Theorem \4-l[ (EW cmd imply that 
the optimal consensus is achieved. 

4.2.2 Proof of Theorem HT^l 

It is easy to find that if 9k = 0, the intersection set in ([1]) is the line segment from Xi{k) to 
Pxi{xi{k)) and then P^"'{k) = Pxi{xi{k)). Then the evolution of the approximate projected 
consensus algorithm becomes 



Xi{k + l)= ^ aij{k)\{l - aj^k)xj{k) + aj^kPxj{xj{k))y i = l,...,n. (31) 
We complete the proof by the following two parts. 
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(i) We first prove that if X^fc^o'^fc ^ then there exist initial conditions from which all 
agents will not converge to Xq. Let z* £ i?"^, which shall be selected later, and Xi{0) = 
z*,i = 1, n. 

By (j3T]) . can be rewritten as 

= Yl ((1 - «i,o)x,(0) + a,, oPx, (xj(0))) 

ieM(o) 

jeM(o) ieM{0) 
= (l-A,o)z* + /3i,oPxo(^*) + Aio, 

where l-^.o = EjGM(o) aij(0)(l-aj,o) and A^o = EjGM(o) aii(0)aj,o(i^x, (^;*) --Pxol^*)) 
with |Ajo| < a^d* for all i. 

We also have 

ieM(i) 

= ^ a,,(l)(l - a,- i)((l - /3,-o)z* + /3,-oPxo(^*)) + 
jeM(i) 

+ Y «ii(l)"i,i^Xo ((1 - /3j,o)2* + Pj,oPxo{z*)^ 
ieM(i) 

= (l-ft,i)z* + A,iPxo(^*) + Aa, 

where 1 — = ^jg_;v'j(i) ~ ~ /^i,o); the third equality follows from Lemma 

EJ (iii) and = A^^V Af^ + Af^ with 

Ail = X] ~ aj,i)Ajo; 

jeM{i) 

Aa= Y «^J(lKl(^'x,(x,(l))-Pxo(a;i(l))); 
jeM(i) 

A?i= 5] a,,{l)aj,i[Pxoixj{l))-Pxo{{l-/3,fi)z* + f3j,oPxo{z*))). (32) 
jeM{i) 

Now we give an estimate of the upper bound of Aji. By Lemma 12.21 (i), 

\Pxo{xj{l)) - Pxoiil - M^* + (3j,oPxo{z*))\ < \xj{l) - {{I - (3j,o)z* + f3j,oPxo{z*))\ 

= lAjol, 

which implies that \Aj^ \ + \A^^\ < maxi<j<„ |Ajo| < a^d* and then |Aji| < \Aj^ \ + |A?;^| + 
lAfJ < {a^ + at)d* for all i. 

Similarly, we can show by induction that for all i and k, Xi{k + 1) can be expressed as 

Xi{k + 1) = (1 - Pi^,,)z* + ft,fcPxo(^*) + A,fc, (33) 



16 



where |Ajfc| < Yl^=Q'^t'^* ^"^^ {Pi,k-,i ^ ''^ifc^o satisfies 

1-A,fc= «*j(fc)(l-ai,fe)(l-/3j,fc-i). (34) 

Based on we can show by induction that 

It follows from ()33p and Lemma 12.21 (ii) that 

|x.(A; + l)Uo > |(l-A,fc)^*+A,fc^'xo(^*)U„-|A.fc|. (36) 

Moreover, for a convex set K and any < A < 1, 

1(1 - A)x + \Pk{x)\^ = 1(1 - A)x + \Pk{x) - Pk{{1 - X)x + XPk{x))\ 
= \{1 - X)x + XPk{x) - Pk{x)\ 

= il-X)\x\K, (37) 

where the second equality follows from Lemma [2^2] (iii). By taking A = /3j_fc in (jST]) . we 
have 

\{l-(3,^k)z* + (^^,kPxo{z*)\x, = (i-A,fc)k*ko- 

Combining the last equality, ([36|) and ([35]) , we obtain 

k k 

\xi{k + l)\x,>\{{l-ai)\z*\x,-Y.aid\ (38) 

/=0 1=0 

Taking the inferior limit on the two sides in (|38p . we have for all 



CXD 



liminf |3;j(A:)|xo > TT(1 - at)\z*\xQ - V'a^d*, 
which is positive provided that 



«=0 Z=0 



l-lxo > i^'7f' (39) 



where n£o(-'- ~ ^t) ^ *-* since < af < 1 for all / and "^^u^o^t ^ Thus, all agents 
can not achieve an optimal consensus for all initial conditions satisfying (j39p . 



(ii) In fact, if X^fclo '^t ^ there is y* = y*{z*) such that limfc_j.oo Xi{k) = y* for all i 
provided z* satisfies ([39|) . 

Denote d = sup;^<j<„ ;j>o \xi{k)\xi- By ([HI]) , we have for all i and k, 

Xi{k + 1)= aij{k)xj{k) +Ti^k, (40) 



17 



where Tj^fc = EjGM(fc) aij{k)ajk{Pxj{xj{k)) - Xj{k)) with \Ti^k\ < da'^ 
Take iq £ V. It follows from (|40p that for all j and A;, 



+ 1) — Xig{k)\ < max |xp(/i:) — Xq{k)\ + max ir^^fcl 



(41) 



Then 



\Xio 

{k + 2) - x,,{k)\ 

ajQj(A; + l)xj(A; + 1) + ^ aioi(^ + l)rj,A:+i — a^jj, 

< max \xj{k + 1) — Xif^{k)\ + max ir^^fc+il 

l<r<n ' 

< max \xp{k) — Xq{k)\ + max iTr^fcl + max Ir^^fc+il 

l^Pi'?^'* l<r<n ' l<r<n ' 

< ^ max |xp(A;) - Xg(A:) I + (i(a^ + a^_^-^). 

We can similarly show by induction that for all k and /, 



k+l-l 

|xj„(/c + /) - Xi(,(A:)| < max \xp{k) - Xq{k)\ + d 

l<p,q<n ^ — ' 

s=k 



< max \xp{k) — Xq{k)\ + dy 



(42) 



Since limfc^oo = and {xj(A:)}^Q is bounded for all i, the consensus is achieved by 
Lemma 14.81 Combining the consensus and the boundedness of {2;j(^)}fc^0' ^^i^re exist 
a subsequence {ki}^Q and y* such that limi^oo Xi{ki) = y* for all i. Therefore, since 
YlT=o'^k ^ ^ > 0) there exists Iq = lo{e) such that \xi{ki) — Xj{ki)\ < e/2 and 

^Yl^ki '^t ^ for all i,j and I > /q- Thus, for / > Iq and /c > ki, 

oo 

|xio(A;) - Xjo(A:/)| < max \xp{ki) - Xq{ki)\ + < e, 

l<p,n<n ^ — ' 



s=k, 



which implies that limi_j.oo Xi^ (k) = y* since e can be arbitrarily small. By what we have 
proven in the first part (i), y* ^ Xq. The desired conclusion follows from the fact that iq 
is taken arbitrarily. 



The proof is completed. 



□ 



5 Critical Angle Error 

The boundedness of system states plays a key role for various optimization methods |28t I26j . 
In this section, we consider the effect of the angle error 9k on the boundedness of the states 
{xi{k),i G V}^o generated by the approximate projected consensus algorithm. 

Suppose ai^k = 1 and 9k = 9 with < 9 < it/2. First the following conclusion shows that 
when 9 < 7r/4, the trajectories of the algorithm are uniformly bounded with respect to all initial 
conditions. 
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Proposition 5.1 Suppose Al, A2, A5 hold and < < 7r/4. Then we have 



suplimsup |3;j(/c)|xo < oo 

a:(0) fc— >oo 



for i = 1, . . . , n. 



Proof: Recall d* = sup^^^ j^^eU" 1'^^ ~ "^^l- We claim that for all i and all initial conditions 
x(0), 

limsup |xi(fc)|xo < 



fc^-oo J- — tany 

which implies the conclusion. 

Based on ([8]) , ([9]) and the definition of d* , we always have 

\Prik)\xo < tane\x,{k)\x, + \PxMimxo 

<tane\xi{k)\xo+d\ (43) 

Furthermore, in the case of \xi{k)\xo > 2d* /{I — tan^), by (03]) we further obtain 

\Pr{k)\xo <tane\xi{k)\xo+d* 

( 1 — tan 0\ , , , , , 

< (l ^ jk.Wko- (44) 

We next consider h(k) := maxi<j<„ \xi{k)\xo by the following two cases: 

(i) h{k) < 2d* /{I - tan 61). By LemmaOand (gS]), we have for ah i, 

\x^{k + l)\xo< E a,,{k)\P^^-mXo 

= ^ Oij (k) [tan 9\x j {k)\ Xo + d*) 
2d* tan d 
2d* 

Therefore, h{k + 1) < 2d* /{I - tanO). 

(ii) h{k) > 2d* /{I - tanO). Define Vi{k) = {i\ \xj{k)\xo < 2d* /{I - ianO) for ah j £ Mi{k)} 
and V2ik) = V\Vi(A;), where V2{k) is nonempty. 

1) If z G Vi{k), by similar analysis with ()45p . we obtain \xi{k + l)|xo ^ 2d* /{I — tanO). 

2) If i G V2ik), we can define MUk) = {j\ j G JVi{k), \xj{k)\x,, < 2d* /{I - ianO)} and 
Mf{k) = J\fi{k)\J\fl (k) , where A/'f(/c) is nonempty. 

Again by Lemma l2.1| we have 

\x,{k + l)\x,< «.i(^)l^'"(^)Uo + J2 a,,{k)\P^'-{k)\x,. (46) 
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For j £ Mlik), by (03]) we have |P/"(A;)|xo < tane\xj{k)\xo + d* < 2d*tan6l/(l - 
tan6l) + d* < 2d* /{I - tan 6*). As a result, by (03]) and (06]) we have 

\xiik + l)|xo < ^ ~ " — ^~~~) ^ aij{k)\xj{k)\xo- 

(47) 

Since \xj{k)\xo > 2d* /{I - tan 6*) for all j G J^f{k), maXjg_/v'i(fc) ki(^)Uo > 2d* /{I - 
tan^). Thus, it follows from ([17|) that 

+ l)|xo < ( V aij{k) + {l-- — V aij{k)) niax |j;i(^)|xo 

/ , , , 1 — tan 6* s \ , , , 
< ((1 - •?) + l{i 2 * + ' 

= (i-5<1^)K*= + i), 

where < 1 — 7/(1 — tan0)/2 < 1 and rj is the lower bound of weights in A2. 
Therefore, based on the two cases (1) and (2), we show that ii h{k) > 2d* / (1 — tan 9) , then 
h{k + 1) < max{2d*/{l - tan 61), (l - rj{l - tan 6)/ 2) h{k)} . 

Combining the two cases (i) and (ii), we have 



h{k + l) < 



2d*/{l-tane), if h{k) < 2d* /{I - tan 9); 

max{2d*/{l - tanO), (l - r/(l - tan 6)/ 2) h{k)} , otherwise 



Thus, the conclusion follows. □ 
We next investigate the non-conservativeness of -7r/4 in Proposition 15.11 We just focus on a 
special case with only one node in the network. Its set is denoted as X^,. We denote the states 
of the node as {x*{k)}'^^Q driven by the approximate projected consensus algorithm: 

x^k + l) £ ^£(x,(A:),e), 

where ^'£{x4k),e) = CxAx*{k),9)f]UxAx*{k)). 
The following conclusion holds. 

Proposition 5.2 Suppose 6 = 7r/4. 

(i) For any compact convex set and any initial condition x*(0) G i?™, we have 

limsup |x*(/c)|x. < \x^{Q)\x,] 

k—^oo 

(ii) There exists an approximate projection sequence {P^°'{k)}'^^Q, P^"'{k) G ^^^(3;^,(A;),7r/4) 
such that 

(11.1) limsup^,^oo |3;*(^)|x, = when X^ is a hall with radius r > 0; 

(11.2) limsupfc_^j^ \x*{k)\x^ = |3;*(0)|x, when X* is a single point. 
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Proof: (i) The conclusion follows from that 

\x^{k + < + 1) — Px,{^*{k))\ < tan9\x^,{k)\x, = \x^,{k)\x,- 

(ii) We just select {Pr(fe)}^o which Ang{P^''{k) - x^{k), Px,{x*{k)) - x,{k)) = tt/A 
for all k satisfying the definition ([1]) such that the solution of the algorithm with ai^k = 1 and 
^fc = 7r/4 satisfies the following two cases. 

(ii.l) It is easy to find that for all k, 

\x,ik + l)\x, = ^J\x.{k)\x.+r^-r, (48) 

which implies that the sequence {|2;*(A;)|x,}fc^o non-increasing and then converges to some 
nonnegative number /x. Letting A; — )• oo in (08]), we have = y^//^ + — r and then = 0. 



x*{k) 



7r/4\ 



r 



/^PxAMk + I)) 



/ r 



Figure 3: Approximate projection with respect to a ball. 

(ii.2) The conclusion is straightforward. □ 
We present another result for the case when 9 > 7r/4 to reveal that, in this case, the node 

states will be unbounded as long as the distance between the initial condition and X* is larger 

than a certain threshold. 

Proposition 5.3 Suppose 9 > 7r/4. Then for any compact convex set X^,, there exists an 
approximate projection sequence {P^"'ik)}^Q such that 

limsup |x*(A;)|x, = oo 
for all initial conditions satisfying |x*(0)|x, > sup^^^ uj2ex, ~ '^2|/(tan0 — 1). 

Proof: Select {Pr(^)}^o ^^^ich Ang(P;"(A:) - x^k), Px, {x^{k)) - x,{k)) = 9 for ah k. For 
all k we have 

\x.{k)\x. > \x4k) - PxAx*{k - 1))| - \Px,{x*{k - 1)) - PxAx*{k))\ 
= tan9\x,{k - l)\x, - \PxAx*{k - 1)) - PxAx*{k))\ 



> tan 6'|x^,(/i: — — sup |a;i — 072! 

UJl,LU2&X, 



(49) 
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At this point, we need the following conclusion: consider a nonnegative sequence {-Zfcj^o 
with Zk^i > (tan0)zfc — d and 6 > 7r/4. Then linifc^oo -^fc = oo if (tan^ — 1)zq — d > 0. Note 
that zi — zq > (tan 9 — 1)zq — d > and then Z2 — zi > (tan 6 — l)zi — d> (tan 9 — 1)zq — d > {). 
We can show similarly by induction that z^+i — > (tan0 — l)zo — for all k. The conclusion 
follows from (|49p and the above conclusion. □ 

Combining Pr op osit ions 1 5 . 1 1 1 5 . 2 1 and 15.31 we see that 7r/4 is a critical value of the angle error 
in the approximate projection regarding maintaining bounded iterative states, li 9 < vr/4, the 
system trajectories are uniformly bounded; if > -7r/4, the trajectories diverge for a special case 
with one single node and particular approximate projection points; \i9 = 7r/4, the trajectories of 
the algorithm with one node are bounded (no longer uniformly with respect to initial conditions) 
and the property of the trajectories highly depend on the shape of the convex set. 

6 Numerical Example 

In this section, we provide a numerical example which shows that the approximate projected 
consensus algorithm presented in this paper may lead to a faster convergence than the projected 
consensus algorithm presented in j?7j . 

Consider a network with three nodes 1, 2 and 3. Suppose Xi{k) £ = 1,2,3. Let Xi, 
X2 and X3 be three unit disks with centers (1,0), (—1,0) and (0,-1), respectively. Set Xi 
corresponds to node i for i = 1,2,3. The communication graph is completely connected. Take 
ail = 0-22 = ^33 = 0.5 and all other weights are 0.25. Here 9k ^ 0, which is a special case of the 
approximate projected consensus algorithm. 

The three agents have the same initial condition 2;i(0) = X2(0) = X3(0). We consider 2500 
initial conditions {(-1.96+0.08p, -1.96+0.08g)| < p < 49, < g < 49} equally spaced over the 
square {(zi, Z2)| |-2i| < 2, |2;2| < 2}. The initial conditions from which the approximate projected 
consensus algorithm with Ui^k = 0.5 converges faster than the projected consensus algorithm 
presented in [27] (compare their distance function h{k) = maxi<j<3 \xi{k)\xo at time 2000) are 
labeled with stars in Figure 4. The decay process of h{k) of a concrete example with initial 
condition a;i(0) = a;2(0) = ^3(0) = (1.8,0.8) is given in Figure 5 (drawn from time 20), from 
which we can see that the approximate projected consensus algorithm has better performance 
from the viewpoint of convergence speed. 

7 Conclusions 

In this paper, we presented an approximate projected consensus algorithm for a network to 
cooperatively compute the intersection of a serial of convex sets, each of which is known only to a 
particular node. We allowed each node to only compute an approximate projection which locates 
in the intersection of the convex cone generated by the current state and all directions with the 
exact projection direction less than some angle and the half-space containing the current state 
with its boundary is a supporting hyperplane to its own set at its exact projection point onto its 
set. Sufficient and/or necessary conditions were obtained for the considered algorithm on how 
much projection accuracy is required to ensure a global consensus within the intersection set. 
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Figure 4: Approximate projected consensus algorithm converges faster than projected consensus 
algorithm for the initial conditions labeled stars. 



0.045 

0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 
0.005 




10 




Figure 5: A typical solution for which approximate projected consensus algorithm (APCA) 
converges faster than projected consensus algorithm (PC A). 
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under the assumption that the communication graph is uniformly jointly strongly connected. We 

also showed that 7r/4 is a threshold for the angle error in the projection approximation to ensure 
a bounded solution for iterative projections. A numerical example was also given indicating that 
the approximate projected consensus algorithm sometimes achieves better performance than the 
exact projected consensus algorithm. This implied that, individual optimum seeking may not 
be so important for optimizing the collective objective. 
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